Localization Transition of Biased Random Walks on Random Networks 
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We study random walks on large random graphs that are biased towards a randomly chosen but 
fixed target node. We show that a critical bias strength b c exists such that most walks find the target 
within a finite time when b > b c . For b < b c , a finite fraction of walks drifts off to infinity before 
hitting the target. The phase transition at b = b c is a critical point in the sense that quantities like 
the return probability P(t) show power laws, but finite size behavior is complex and does not obey 
the usual finite size scaling ansatz. By extending rigorous results for biased walks on Galton- Watson 
trees, we give the exact analytical value for b c and verify it by large scale simulations. 
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Random walks are a fascinating subject, both for their 
intrinsic mathematical beauty and for their wide range 
of applications One of the most celebrated results 
is that unbiased random walks on regular lattices are 
recurrent for dimension d < 2, while they are transient 
for d > 2 0. This means that a walk starting at a node 
x certainly returns to x when d < 2, but has a finite 
chance to escape to infinity when d > 2. 

Here we study an analogous problem for walks on ran- 
dom graphs which lack small loops in the limit of infinite 
graph size. Examples of such graphs include the Erdos- 
Renyi (ER) random graphs, as well as random graphs 
with any fixed degree sequence, provided that the vari- 
ance of the degree distribution is finite. We call the lat- 
ter Molloy-Reed (MR) graphs 0. An ER graph with N 
nodes is constructed by introducing a link between each 
pair of nodes with probability p. It is in many ways sim- 
ilar to an infinite dimensional lattice. In particular, its 
diameter I increases only logarithmically with the num- 
ber of nodes, while N ~ l d on a d-dimensional lattice. 
Thus one should expect unbiased random walks on an 
ER graph to be transient. Although less is known rigor- 
ously about MR graphs, we expect the same to be true 
for them. In order to arrive at a non-trivial problem, 
we thus consider walks biased towards a randomly cho- 
sen but fixed "target" node. We show that there is a 
phase transition from recurrence (or localization) to de- 
localization at a critical bias strength. Notice that this is 
unrelated to similar phase transitions observed e.g. in 
where the bias is not towards but away from the target. 

This seemingly abstract mathematical problem has a 
number of practical applications. Consider, e.g., rout- 
ing a message from node A to node B on the internet. 
Each node i maintains a routing table which indicates for 
each target node the optimal first step, parting from i. 
Using her routing table, A sends the message to her opti- 
mal neighbor i\. From there it is sent, using the routing 
table at i\, to i%, etc., until i n = B is reached. If all 
routing tables are correct and up-to-date, the message 
reaches its destination along the optimal path (i.e., the 
path between A and B with the shortest length). How- 
ever, some routing tables may be faulty, either because 



they contain mistakes or have become obsolete due to 
changes in the internet topology. If the fraction of such 
nodes is below a certain threshold, the effect is small and 
routing is still efficient: the average time to go from A 
to B scales linearly with the distance. But when the 
fraction exceeds a critical value, the message might take 
a very long and convoluted path through a finite frac- 
tion of the entire network, before reaching its destina- 
tion. There are of course many details in which routing 
on the real internet differs from the simpler problem of 
biased random walks on random graphs [5|. For instance, 
the degree distribution of the internet is approximately 
scale- free with a divergent second moment [y] ; stray mes- 
sages on the internet are killed after some time; and the 
bias is implemented differently and is quenched (the same 
routing tables are used at successive time steps). Despite 
these technical details, the two problems are basically the 
same. For a related discussion of communication based 
on noisy routing, see 

Another application is to quantum mechanics. A ran- 
dom walker on a graph corresponds quantum mechani- 
cally to a particle with hopping dynamics. The recur- 
rence problem for an unbiased walk corresponds then to 
the question whether such a particle, subjected to an at- 
tractive (5-potential concentrated on a single node, forms 
a bound state. Localization of biased random walks cor- 
responds to the existence of bound states for potentials 
which increase linearly with distance from the target 
node. Derealization would imply the paradoxical sit- 
uation that no bound state exists, although the potential 
increases forever as one goes further away from the tar- 
get node. Instead, a particle released near the bottom of 
the potential continues to climb up the potential forever, 
because there are always more paths leading uphill than 
leading back to the bottom. 

Technically, we define our walks as follows. Consider a 
finite but large undirected random graph Q with N nodes 
and with degrees chosen from a distribution Pg(k). We 
assume that Pg(k) is such that Q is sparse and has a 
giant connected component (GC). We define the distance 
between any two nodes as the number of links in the 
shortest path connecting them. We randomly choose a 
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node A on Q and label all other nodes according to their 
distance from A. Consider node i at distance di from A, 
which has k~ , and neighbors at distances di + 1, 
di — 1, and dj from A, respectively. The random walk 
steps from i with probabilities 
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to a node closer to, at the same distance as or further 
away from A, respectively. Finally, normalization re- 
quires Mi = bk^ + kf + b~ x kf . 

Similar walks on regular lattices have been studied re- 
peatedly, see e.g. ||. More importantly, Eq. ([1]) is a 
generalization of the bias used in the 'A-biased random 
walks' studied on Galton- Watson (GW) trees SfMEII- 
Starting at the root, each node on the tree has a number 
of daughter nodes chosen from a prescribed distribution. 
The root is chosen as the target A of the random walk. 
On a tree, every node i ^ A has only one neighbor closer 
to the target, k~ = 1, and no neighbors at the same 
distance, kf = 0. The probabilities of the next step are 
chosen such that 0, [lj| Pi IpI — A- This corresponds to 
Eq. JT]), restricted to tree-like graphs, with b = \f\. 

The typical length of a loop in an undirected graph 
Q with finite mean degree and finite degree variance is 
of order In TV [H, EH- If only local properties are of 
interest, the graph can, as N — > oo, be effectively re- 
placed by a GW tree. However, a subtlety should not 
be overlooked: to obtain a rooted tree T from a loopless 
undirected graph Q, we have to choose randomly a node 
A as the root of the tree, and draw an arrow on each 
link pointing away from A. A node with degree k on Q 
has in-degree 1 and out-degree k — 1 on the effective tree. 
The out-degree distribution on the tree can be related to 
the degree distribution Pg(k) of the graph Q, 
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with k = k + + 1. The prefactor k/(k) on the right hand 
side takes into account the fact that each of the k links 
attached to the node can play, with the same probability, 
the role of the incoming link on the tree. 

A GW tree will grow to infinity only if its average 
out-degree is larger than 1 [14]. This implies that the 
underlying random graph is connected when 0, [l5| 



M= (fc + )GW = (Hk-l))g/(k)g > 1, 
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which we assume to be satisfied. Further for the graph to 
be sparse and lack small loops we assume that /x <C N. 

If there is a critical bias A c for walks on a GW tree 
so that walks are localized near the root for A > A c and 
delocalized for A < A c , then an analogous critical bias 
b c < \/Ac must also exist for graphs. The reason is simply 
that only local properties arc relevant when the walks are 
localized [l8[. Hence localization on a GW tree implies 
localization on the corresponding graph. Conversely, as- 



sume that b c were strictly smaller than ^f\~ c . Then walks 
on the GW tree with b 2 < A < A c would also be lo- 
calized, contradicting the starting assumption that the 
critical point is at A c . Thus we must have b c = \/A7. 

Indeed, it is known that the A-biased random walk on 
a GW tree has a localization transition at A c = /i, such 
that the walk is recurrent for A > A c and transient for 
A < A c 11]. For A = A c the authors of 0] prove a central 
limit theorem which states that the walk behaves, as far 
as the distance from the root is concerned, like unbiased 
1-d Brownian motion. In the following we will summarize 
the arguments leading to those conclusions and indicate 
modifications required to apply to ER and MR graphs. 

For a walker at a node i that is different from the root 
A, we define pf to be the probability that it escapes to 
infinity before it hits the root. We denote by V{i) the 
set of parents of i, i.e. the neighbors that are closer to 
A than i. Similarly, S(i) is the set of siblings of i (i.e. 
dj = di for all j S S(i)), and C(i) denotes its children. 
Finally, we define = 0. Then we have, for any i ^ A, 
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Using Mi defined after Eq. ((T|) we can rearrange this to 

jec(i) 



jer(i) ' jec{i) ' jes(i) 

(5) 

A sum over all i with fixed di = d, cancels all the contri- 
butions of the siblings. Defining 



we get the recursion 
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which can be iterated to give the average escape proba- 
bility from all children of the root 
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For trees this simplifies because each node (except for the 
root) has only one parent, and the number of terms on 
the r.h.s. increases in average as fi d for large d. Since 
each term is bounded, the total sum increases at most as 
([i/b 2 ) d . Thus pe = for b 2 > //, showing that the walk 
is recurrent. The proof that b c is not only < ^JJi, but 



b c = yfji , 



(9) 



is found in [llj. Basically, one replaces in Eq. ([5]) the 
number of terms by its expected value and the difference 



3 



pf - pf by its average (p)d+i - (p)d, to obtain 



Pe ~ (^r) ((p)d+i - (p)d) 



which has the solution 



(p) d =a- /3 



(10) 
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with a and (3 being constants. For b 2 > p, the only 
consistent values are a — /3 — 0. For b 2 < p one has a 
non-zero solution, indicating p^ > 0. 

For general finite MR graphs the number of nodes at 
distance d from the root increases slower than fj, d , due to 
the existence of loops. In this case b c < ^fjl holds a for- 
tiori, but the rigorous proof that also b c > JJi becomes 
more difficult. The corrections to \i d are small for small 
d (there are few small loops) , and we can expect that the 
effect of loops on b c can be neglected in the limit N —> oo. 
But loops should effect the finite size scaling behavior. 

To verify the existence of a localization transition on 
large ER graphs, check Eq. ([9]), and study finite size ef- 
fects, we perform large scale numerical simulations. We 
first generate large ER graphs (with up to 2 x 10 8 nodes), 
for several nominal values of (k). These are chosen such 
that the GC contains more than 90% of all nodes. Prun- 
ing all nodes and links not connected to the GC gives 
the final graph size N. This procedure increases (k) 
slightly and makes the degree distribution slightly non- 
Poissonian, so that fx no longer coincides with the naive 
theoretical estimate fi = (k) for ER graphs. Instead, it 
is estimated from the exact definition Eq. ([3]). 

After extracting the GC, we choose randomly one of 
its nodes as the target A and compute every other node's 
distance from A. Then we start M walks at A and follow 
them until they return to A. This is repeated M' times by 
taking new target nodes, and finally the whole procedure 
is repeated for M" different graphs. 

The first observable to be discussed is the average re- 
turn time. For b = 1, i.e. for unbiased random walks, the 
average return time (t) on a finite graph with TV nodes is 
oc N [12], We expect the same to hold for all 1 < b < b c . 
For the critical bias b — b c , since the distance of the 
biased walk from the root behaves as a 1-d Brownian 
motion, the return time should grow as logiV which is 
the diameter of the graph. In contrast, linijv— >ao(i) < oo 
in the localized regime b > b c . 

In Fig. [1] we show (t) versus N for various biases on 
a linear-log scale. The nominal number L/N of links 
per node in constructing the ER graphs was 7/3, which 
would give ,u = 14/3 = 4.6667. For the GC, we found 
numerically fi = 4.6673. The critical bias should be b c = 
2.1604, which certainly agrees with with Fig.[TJ 

More significant is the distribution of return times. In 
the inset of Fig. [5] we plot (not normalized) histograms of 
return times P(t) against t for walks with various values 
of b, on the GC of an ER graph with N = 2 x 10 8 and 
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FIG. 1: (color online) Average return times against the sys- 
tem size N for different values of the bias b, on ER graphs 
with (k) = 14/3. 
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FIG. 2: (color online) Histograms of return times for walks 
on ER graphs with N = 2 x 10 8 and L/N = 1.5. Shown is 
t 3 ^ 2 P(t) where P(t) is the return probability. The inset is the 
plot of P(t). The critical curve (third curve from above) is 
the one which is flattest at t ~ 100. Each curve corresponds 
to a fixed bias b, with steps of 0.04. The upper (red) curves 
correspond to even t, the lower (blue) ones to odd t. 



L/N = 1.5. For this graph b c = 1.73542. The most strik- 
ing observation is a very large even/odd effect for small 
t: Compared to the values for even t, Pit) is nearly zero 
for odd t, an immediate consequence of the suppression 
of small loops. On a loopless tree, return times are al- 
ways even. Conversely, we can assume that loops are 
negligible for even times for which P(t ± 1) <C P(t). In 
this regime we should thus expect that the result for GW 
trees applies, i.e. P(t) ~ t~ 3 / 2 at the critical bias. This 
prediction (which is the same as for unbiased 1-d ran- 
dom walks |9]) is in complete agreement with our data, 
as seen from Fig. [21 Fig. [2] also indicates the complicated 
finite-size behavior. P(t) is not convex at the critical 
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FIG. 3: (color online) Density of walkers against distance 
from the target, for five values of b close to criticality. 



point, and the usual finite size scaling ansatz (power law 
times homogeneous function) does not hold. 

Finally, if critical walks resemble unbiased 1-d Brow- 
nian motion, with reflecting boundaries at d = and at 
d w IniV 0], then the density p(d) of walkers per launch, 
integrated over all times, should be constant at b — b c . 
Fig. [3] shows that p{d) is indeed flat, to very high accu- 
racy, for b = 1.732 ±0.002. This agrees with the exact b c 
within two standard deviations, and gives the most pre- 
cise numerical verification of the theoretical prediction. 

We have shown that walks on a random graph biased 
towards a given node show a localization/delocalization 
transition, corresponding to a transition from recurrent 
to transient behavior. If routing is noisy but the noise 
level is below a critical threshold, the walks are able 
to reach their destination in a finite amount of time. 
This resilience of the walks can be taken into account 
while designing routing strategies for the internet or 



other traffic problems with noisy dynamics. The localiza- 
tion/transition should also be related to a curious effect 
in quantum mechanics on random graphs with linearly 
rising potentials, where we expect a paradoxical "unbind- 
ing" transition when the potential becomes too shallow. 

Our numerical studies have only included ER graphs, 
mainly because these can be easily generated even with 
very large sizes. However, we expect the existence of 
a derealization transition to be robust and to hold for 
any degree distribution with a finite second moment. For 
scale free graphs, P(k) ~ 1/fc 7 with 7 < 3, the second 
moment diverges with N, leading formally to b c = 00. 
For scale free graphs, the number of nodes within dis- 
tance d from the target does not grow exponentially, 
~ yu d as assumed in Eq. (fit)]) . It rather grows super- 
exponentially as exp(exp(d)) [lj| [T3|- To compensate 
this, the critical bias has to be distance-dependent, grow- 
ing also super-exponentially with d. Alternatively, if we 
describe the finite-size behavior by an ./V-dependent but 
d-independent effective critical bias b c (N), then this has 
to grow as a power of N. 

Finally, we should point out that our theoretical pre- 
diction for b c holds only for a very particular type of bias, 
given by Eq. JT]). For other biases we expect the transi- 
tion to show the same scaling laws (as long as the bias 
strength is independent of d), although we can no longer 
predict the exact location of the phase transition. The 
same is true for graphs with many small loops, including 
those often observed in real world networks. If the bias 
strength increases (decreases) with d beyond limit, then 
one has always (de)-localisation for graphs with a finite 
second-moment of the degree distribution and with expo- 
nential growth of the number of neighbors with distance. 
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